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Abstract —We consider the problem of voltage stability and 
reactive power balancing in islanded small-scale electrical net¬ 
works outfitted with DC/AC inverters (“microgrids”)* A droop¬ 
like voltage feedback controller is proposed which is quadratic 
in the local voltage magnitude, allowing for the application of 
circuit-theoretic analysis techniques to the closed-loop system. 
The operating points of the closed-loop microgrid are in exact 
correspondence with the solutions of a reduced power fiow 
equation, and we provide explicit solutions and small-signal 
stability analyses under several static and dynamic load models. 
Controller optimality is characterized as follows: we show a one- 
to-one correspondence between the high-voltage equilibrium of 
the microgrid under quadratic droop control, and the solution 
of an optimization problem which minimizes a trade-off between 
reactive power dissipation and voltage deviations. Power sharing 
performance of the controller is characterized as a function of 
the controller gains, network topology, and parameters. Perhaps 
surprisingly, proportional sharing of the total load between 
inverters is achieved in the low-gain limit, independent of the 
circuit topology or reactances. All results hold for arbitrary 
grid topologies, with arbitrary numbers of inverters and loads. 
Numerical results confirm the robustness of the controller to 
unmodeled dynamics. 

Index Terms —Inverter control, microgrids, voltage control, 
Kron reduction, nonlinear circuits 

I. Introduction 

The wide-spread integration of low-voltage small-scale re¬ 
newable energy sources requires that the present centralized 
electric power transmission paradigm to evolve towards a more 
distributed future. As a flexible bridge between distributed 
generators and larger distribution grids, micro grids continue to 
attract attention 0-0 Microgrids are low-voltage electrical 
distribution networks, heterogeneously composed of distributed 
generation, storage, load, and managed autonomously from the 
larger primary grid. While often connected to the larger grid 
through a “point of common coupling”, microgrids are also 
able to “island” themselves and operate independently 0, 
Q. This independent self-sufficiency is crucial for reliable 
power delivery in remote communities, in military outposts, 
in developing nations lacking large-scale infrastructure, and 
in backup systems for critical loads such as hospitals and 
campuses. Energy generation within a microgrid can be quite 
heterogeneous, including photovoltaic, wind, micro-turbines, 
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Fig. 1. (a) Schematic of a “parallel” microgrid, in which several inverters 

supply power to a distribution bus (effectively a single load) (b) A simple 
non-parallel microgrid consisting of flve loads • and three inverters ■. 


etc. Many of these sources generate either variable frequency 
AC power or DC power, and are interfaced with a synchronous 
AC grid via power electronic DC/AC inverters. In islanded 
operation, at least some of these inverters must operate as grid¬ 
forming devices. That is, control actions must be taken through 
them to ensure synchronization, voltage stability, and load 
sharing in the network 0>0.0>0. and to establish higher- 
level objectives such as frequency regulation and economic 
dispatch 

A. Literature Review 

The so-called droop controllers (and their many derivatives) 
have been used with some success to achieve primary control 
goals such as stability and load sharing, see 0 ^ & 0 - 

GD Despite being the foundational technique for networked 
operation of islanded microgrids (Figure [^, the stability and 
fundamental limitations of droop-controlled microgrids has 
only recently begun to be investigated from a rigorous system- 
theoretic point of view 0 GD-Gg 

Our focus here is on voltage control, which we now provide 
some context for. In high-voltage networks, the grid-side volt¬ 
ages of transformers which interface synchronous generators 
are regulated to nominal set points via automatic voltage 
regulator systems. However, if this strategy was applied to 
inverters in low-voltage networks, the small mismatches in 
voltage set points would combine with the large impedances 
presented by distribution lines, and result in large circulating 
reactive currents between inverters. Moreover, due to the close 
electrical proximity of devices in small microgrids, and the 
small power capacities (ratings) of the sources, it is desirable 
that controllers establish set percentages of the total load to 
be supplied by each inverter (so-called power sharing). These 
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technical obstacles motivated the use of voltage droop control 
fT^ as a heuristic proportional controller to establish power 
sharing between units while maintaining voltages within a 
reasonable range around their set point values; a technical 
review of the voltage droop controller is presented in Section 

IIEb] 

The widespread deployment of voltage droop has led to 
several attempts at stability analysis (Tg,(Tg. Both references 
begin with all-to-all Kron-reduced network models which do 
not explicitly contain loads. They then assume the existence of 
a system equilibrium, and derive sufficient stability conditions 
which depend on the uncharacterized equilibrium. Results 
apply only for constant-impedance or constant-current load 
models, and no characterizations are given of reactive power 
sharing in steady state. Thus, the available literature offers 
no guidance on the foundational issue of operating point 
feasibility, i.e., the existence and locations of steady states for 
the microgrid which satisfy operational constraints. This gap 
of knowledge means that precise reactive loading limits and 
security margins are unknown, making system monitoring and 
non-conservative operation difficult. Moreover, stable network 
operation is often limited by stiff constant power loads, which 
do not reduce their current consumption when voltages fall. 

Regarding control performance, a well-recognized drawback 
of voltage droop control is that — while achieving better 
power sharing performance than mere voltage regulation — the 
power sharing properties can still be quite poor. This has led 
some authors to investigate alternative strategies to implement 
reactive power sharing |[7|, | [TT| , fTT) , |T^ that are based on 
directly load measurements or inter-unit communication. Aside 
from these additional requirements, these approaches do not 
assist in characterizing the power sharing properties of the 
standard, decentralized droop control. A general discussion 
of fundamental limits for voltage control and reactive power 
sharing can be found in Section III], and some recent non¬ 
droop-based control (but nonetheless related) approaches to 
reactive power compensation can be found in p^-|[28|. 

In summary, two fundamental outstanding problems regard¬ 
ing the voltage stability of droop-controlled microgrids are (i) 
to establish conditions under which a stable network equilib¬ 
rium exists for various load models, and (ii) to characterize the 
resulting power sharing properties in terms of grid topology, 
branch admittances, loads, and controller gains. 


B. Contributions 

In the preliminary version of this work p^ , we introduced 
the quadratic droop controller, a slight modification of the 
conventional linear droop curve to a quadratic one. Quadratic 
droop leads to the same steady states as conventional droop 
control (for particular gains) but allows for circuit-theoretic 
techniques to be leveraged for analysis purposes. In p9| we 
studied system stability for the case of an islanded parallel 
microgrid, where several inverters feed power to a single 
common load or a common distribution bus (Remark]^. In this 
paper we depart from the case of a parallel microgrid and study 
system stability for arbitrary interconnections of inverters and 
loads, under various load models. In addition, we investigate 


the optimality and performance characteristics of the controller. 
Linearization ideas presented in pQ| are also relevant to the 
current presentation at a technical level, and will be used in 
proofs of several results. 

There are three main technical contributions in this paper. 
First, in Section HI we present the quadratic droop controller, 
highlighting its circuit-theoretic foundations. We state and 
prove a correspondence between the equilibrium points of 
the closed-loop system and the solutions of a reduced power 
flow equation, before proceeding to analyze the stability of the 
closed-loop system for both static and dynamic load models. In 
particular, we analyze system stability for the static ZI and ZIP 
load models, and for a dynamic model describing a variable 
shunt susceptance. Our stability analyses in Theorem 
Theorem [33) and Theorem [T4| result in succinct and physically 


intuitive stability conditions, significantly generalizing the 
results and intuition developed for parallel microgrids in p9| . 

Second, in Section |IV-A| we investigate the optimality 
properties of the quadratic droop controller. We demonstrate 
that the microgrid dynamics under quadratic droop control can 
be interpreted as a decentralized primal algorithm for solving 
a particular optimization problem, leading to an additional 
interpretation of the quadratic droop controller as an algorithm 
which attempts to minimize a trade-off between total reactive 
power dissipation and voltage deviations. We comment on 
possible extensions and generalizations of these results. 

Third and finally, in Section |IV-B| we investigate the power 
sharing properties of our controller, and closely examine two 
asymptotic limits. In the limit of high controller gains, we 
show that the reactive power is supplied by an inverter to a 
load in inverse proportion to the electrical distance between 
the two, independent of relative controller gains. Conversely, in 
the low-gain limit we find that inverters provide reactive power 
in proportion to their controller gains, independent of electrical 
distance. When neither limit is appropriate, we provide a 
general formula and quantify the power sharing error as a 
function of the controller gains and network parameters. These 
predictions along with the robustness of the approach are 
studied numerically in Section [V] Remarks throughout the 
paper provide details on extensions and interpretations of our 
results, and can be skipped on a first reading without loss of 
continuity. 

At a more conceptual level, this work has two major 
features which distinguish it from the primary literature on 
droop control. The first is our analytical approach, based on 
circuit theory |3T) and approximations of nonlinear circuit 
equations p0| . While we make specific modeling assumptions, 
they are standard in the field 0 and allow us to derive 
strong yet intuitive results regarding stability and optimality of 
operating points. Our focus is on characterizing the operating 
points arising from decentralized voltage control, along with 
characterizing controller performance at the network level. We 
make use of different load models in different sections of the 
paper, in part to simplify certain presentations, and in part to 
illustrate a variety of analysis techniques which should prove 
useful for other researchers. The second distinguishing feature 
is our explicit consideration of loads which are not collocated 
with inverters. Voltages at non-collocated loads are typically 
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lower then the controlled voltages at inverters, and hence it is 
these load voltages which ultimately limit network stability. 

The remainder of this section recalls some basic notation. 
Section |IJ reviews basic models for microgrids, loads, and 
inverters, along with the conventional droop controller and the 
relevant literature. Our main results are housed in Sections 


future research. 

C. Preliminaries and Notation 

Sets, vectors and functions: We let M (resp. M>o) denote the 
set of real (resp. strictly positive real) numbers, and let j be the 
imaginary unit. Given x [x] G is the associated 

diagonal matrix with x on the diagonal. Throughout, and 
are the n-dimensional vectors of unit and zero entries, and 
© is a matrix of all zeros of appropriate dimensions. The nxn 
identity matrix is A matrix W e is row-stochastic 

if its elements are nonnegative and Wlm = In- A matrix 
^ ^ I^nxn is an M-matrix if Aij < 0 for all i j and all 
eigenvalues of A have positive real parts. In this case —DA 
is Hurwitz for any positive definite diagonal D G (so- 

called T)-stability), and A~^ > © component-wise, with strict 
inequality if A is irreducible p^ . 

II. Modeling and Problem Setup 

We briefly review our microgrid model and recall the con¬ 
ventional voltage droop controller. 

A. Review of Microgrids and AC Circuits 

a) Network Modeling: We adopt the standard model of a 
quasi-synchronous microgrid as a linear circuit represented by a 
connected and weighted graph G(V, f), where V = {1,..., n-f- 
m} is the set of vertices (buses) and f C V x V is the set 
of edges (branches). There are two types of buses: n > 1 
load buses C = n}, and m > 1 inverter buses I = 

{n -f 1,...,n + m}, such that V = CUX. These two sets of 
buses will receive distinct modeling treatments in what follows. 

The edge weights of the graph are the associated branch 
admittances i/ij = gij ibij G C[^ The network is con¬ 
cisely represented by the symmetric bus admittance matrix 

Y G ^ where the off-diagonal elements are 

Yij = Yji = —yij for each branch {i, j} G £ (zero if {i^j} ^ 
£), and the diagonal elements are given by Yu = Vij- 

To each bus i G V we associate a phasor voltage Vi = 

and a complex power injection Si = Pi Y}Qi. For dominantly 
inductive lines, transfer conductances may be neglected and 

Y = ]B is purely imaginary, where B G is 

the susceptance matrix. We refer to Remark for a discussion 
of this assumption. The power flow equations then relate the 
bus electrical power injections to the bus voltages via 

Pe,i = BijEiEj sin(6>i - Oj), i G V , (la) 

E Ti-\-Tn 

BijEiEj cos(6>i - Oj) , i G V. (lb) 


III-V] In Section VI we summarize and provide directions for 


In this work we focus on dynamics associated with the reactive 
power flow equation and refer the reader to 0-® (D 
and the references therein for analyses of microgrid active 
power/frequency dynamics. We will work under the standard 
decoupling approximation, where \0i — 0j\ ^ 0 and hence 
cos(6>i — Oj) ^ 1 for each {i, j} G £\ see |[^, |[^. This can 
be relaxed to non-zero but constant power angles Oi — Oj = 
lij < 7r/2 at the cost of more complicated formulae, but we 
do not pursue this here; see Remark Under the decoupling 
assumption, the reactive power injection QQ^i{0,E) in •dUi 
becomes a function of only the voltages E, yielding 

Qe,i{E) = ~Ei ^ EijEj , ( 2 ) 

or compactly in vector notation as 

Q = -[E]BE, (3) 

where E = (F^i,..., For later reference, some well- 

known properties of the susceptance matrix are recorded in 
Lemma [a. 11 

b) Load Modeling: When not specified otherwise, we 
assume a static load model of the form Qi{Ei), where the 
reactive power injection Qi is a smooth function of the supplied 
voltage Ei. With our sign convention, Qi{Ei) < 0 corresponds 
to an inductive load which consumes reactive power. Specific 
static and dynamic load models will be used throughout the 
paper, and we introduce these models when needed. 

c) Inverter Modeling: A standard grid-side model for a 
smart inverter is as a controllable voltage source 

TiEi = Ui, iel, ( 4 ) 

where Ui is a control input to be designed, and > 0 is a 
time constant accounting for sensing, processing, and actuation 
delays. The model assumes that the control loops which 
regulate the inverter’s internal voltages and currents are stable, 
and that these internal loops are fast compared to the grid-side 
time scales over which loads change. This model is widely 
adopted among experimentalists in the microgrid field 
p^ , p7| , and further explanations can be found in 0 and 
references therein. We note that grid-feeding (i.e., maximum 
power point tracking) photovoltaic inverters or back-to-back 
converters interfacing wind turbines can be modeled from the 
grid side as constant power sources, and are therefore included 
in our framework as loads. 

B. Review of Conventional Droop Control 

The voltage droop controller is a decentralized controller for 
primary voltage control in islanded microgrids. The controller 
is a heuristic based on the previously discussed decoupling 
assumption, and has an extensive history of use; see 0,0, 
(n)> (n), oz). For the case of inductive lines, the droop 
controller specifies the input signal Ui in 0 as the local 
feedback Chapter 19Q 

Ui =-{Ei - E*) - niQe,i{E), iel, (5) 


^For the short transmission lines of microgrids, line charging and leakage 
currents are neglected and branches are modeled as series impedances 
Shunt capacitors/reactors will be included later as load models. 


*For grid-connected inverters, one sometimes sees this formula augmented 
with power set points Ui = —{Ei — E*) — ni{Qe^i{E) — Here we 

consider islanded operation where = 0. 
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where £’* > 0 is the nominal voltage for the inverter, and 
Qi{E) is the measured reactive power injection; see 0^(D 
for details regarding measurement of active/reactive powers. 
The controller gain > 0 is referred to as the droop coefficient. 
From ([^, it is clear that if the inverter reactive power injection 
(5e,i(^) is non-zero, the voltage Ei will deviate from E^. 

Remark 1: (Comments on Modeling, Decoupling, and 
Droop Control) Formally, the decoupling assumption leading 
to ^ is an assumption about the sensitivity of reactive power 
injections (3e,i with respect to phase angles Oi — Oj. From 
this sensitivity is proportional to sm{0i —Oj), and around 
normal operating points is therefore roughly zero; phase 
angles effect reactive power only though second-order effects. 
The strong relationship between reactive power and voltage 
magnitudes is why the former is often used as a tool to 
regulate the latter. The voltage droop controller 0 is designed 
under the assumption that the susceptance to conductance 
ratios bij /gij are large within the microgrid. This assumption 
is typically typically justified in engineered settings, as the 
inverter output impedances are controlled to dominate over 
network impedances giving the network a strongly inductive 
characteristic p9| , Chapter 7]. For dominantly resistive (or 
capacitive) microgrids, the appropriate droop controllers take 
different forms Chapter 19], with the following general 
structure: assuming all lines have uniform reactance/resistance 
ratios bijjgij = const., the controller is modified by replacing 
Qe,i in 0 with Qe,i = cos((/9)Pe,i + si^(/:?)(5e,i, whcrc 
tan((/9) = —bij I gij. A standard calculation yields (5e,i = 

— Y^^=r^EijEiEjCos{Oi — Oj), where B is an admittance 

matrix with edge weights —{b‘fj This formula is 

of the same mathematical form as ( |lb| ), and for such uniform 
networks, we can restrict ourselves — without loss of generality 

— to the decoupled reactive power flow 0 and the conventional 

droop controller 0. As a special case of the above, we note 
that all analysis herein also applies to voltage/active power 
droop control for dominantly resistive microgrids. □ 

Despite its extensive history, the voltage droop controller 
0 has so far resisted any rigorous stability analysis. In our 
opinion, the key obstacle has been — and remains — the 
difficulty in determining the high-voltage equilibrium of the 
closed-loop system for arbitrary interconnections of devices. 
We remove this obstacle in the next section by proposing a 
controller based on the nonlinear physics of AC power flow. 


III. Quadratic Droop Control 


A. Definition and Interpretation 

While the droop controller ^ is simple and intuitive, it is 
based on the linearized behavior of AC power flow around 
the system’s open-circuit operating point, and does not respect 
the inherently quadratic nature of reactive power flow. We 
instead propose a physically-motivated modification of the 
conventional voltage-droop controller In place of 0> 
consider instead the quadratic droop controller 


Ui = KiEi{Ei - E;) - Q,4E ), i G X, 


( 6 ) 


where Ki < 0 is the controller gainj^ Compared to the conven¬ 
tional controller 0, the gain on the regulating term {Ei — E’*) 
in ([^ now scales with the inverter voltage. Combining the 
power flow 0 with a static load model Qi{Ei) at each load 
bus i e jC, we must also satisfy the power balance equations 


0 = Qi{Ei) - Q,^i{E), ieC. (7) 


Combining now the inverter model 0, the quadratic droop 
controller 0> the load power balance 0, and the power flow 
equation ([^, the closed-loop dynamical system is differential- 
algebraic, and can be written compactly as 



f Ql{El) a 
\[Ej]Kj{Ej-E})) 


E[E]BE, 


( 8 ) 


where Ej — (F/^^-i-i, • • •, Ql{El) — 

(Qi(X'i),..., Qn{En)), and r/ and Kj are diagonal matrices 
with elements and Ki respectively. Since the variables Ei 
represent voltage magnitudes referenced to ground, they are 
intrinsically positive, and for physical consistency we restrict 
our attention to positive voltage magnitudes. 

Remark 2: (Interpretations) Before proceeding to our 
main analysis, we offer two interpretations of the quadratic 
controller 0. one theoretical and one pragmatic. 

Circuit-theoretic interpretation: The design ([^ can be inter¬ 
preted as control by interconnection, where we interconnect the 
physical electrical network with fictitious “controller circuits” 
at the inverter buses i G X pO] . Indeed, consider the two-bus 


E^ 


K, 


Fig. 2. Linear circuit representation of the quadratic droop controller 0. 


circuit of Figure where the blue bus has variable voltage 
Ei and is connected via a susceptance < 0 to the green 
bus of fixed voltage X’*. The current-voltage relations for this 
fictitious circuit are 



where li (resp. /*) is the current injection at the bus with 
voltage Ei (resp. X*). Now, let there be m of the two-bus 
circuits in Figure one for each inverter. If we identify the 
variable-voltage blue buses of these circuits with the inverter 
buses of our original network, and impose that the current 
injected into the former must exit from the latter, we obtain an 
augmented network with n + 2m buses, and in vector notation 
the current-voltage relations in the new network are 

(Il\ (Bll Bli © \ /EA 

I ©m = J I Bji A Ki —Kj I X/ . (10) 

\rjJ V© -Ki Ki J \EiJ 

where we have block partitioned all variables according to 
loads, inverters, and fictitious controller buses. In this aug¬ 
mented circuit, the inverters behave as interior nodes joining 
the fictitious controller buses to the loads, and do not sink or 

^ While the chosen sign convention Ki < 0 is unconventional, it will 
simplify the resulting formulae and help us interpret the controller physically. 
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Fig. 3. Diagram showing network augmentation and reduction. First, each inverter bus a of the network in F igure |l(b)| is interconnected with a two-bus 
controller circuit, consisting of an inverter bus B and fictitious bus ♦ at constant voltage . In Theorem |3.1| the inverter buses are eliminated via Kron 
reduction, leaving a reduced network with only fixed voltage buses ♦ and load buses •. 


source power themselves. Left multiplying the first two blocks 
of equations in ( p^ by [E] and noting that by the definition 
of reactive power ^El]Il = Ql{El), we immediately obtain 
the right hand side of Thus, the closed-loop equilibrium 
points under quadratic droop control can be interpreted as the 
solutions of the power balance equations for an expanded linear 
circuit. 

Practical interpretation: Far from being linear, volt¬ 
age/reactive power capability characteristics of synchronous 
generators display significant nonlinearities. In the absence of 
saturation constraints, the characteristics are in fact quadratic 

f Equation 3.105], and thus the quadratic droop controller 
more accurately mimics the behavior of a synchronous 
generator with automatic voltage regulation, compared to the 
classical controller 0- This quadratic dependence means that 
the marginal voltage drop (voltage drop per unit increase in 
reactive power) increases with reactive power provided. □ 
Remark 3: (Generalizations) The quadratic droop con¬ 
troller ([^ is a special case of the general feedback controller 

Ei = Ei (oLijEj -h PijEj^ , (11) 

where aij = aji and /3ij are gains. One recovers 0 by setting 
an = Ki, I3ii = —Ki, and all other parameters to zero. While 
the decentralized controller ([^ can be interpreted as control-by- 
interconnection with the circuit of Figure the controller ( pTj ) 
represents a more general, densely interconnected circuit with 
m variable-voltage nodes and m fixed voltage nodes. Since 
decentralized control strategies are preferable in microgrids, we 
focus on the decentralized controller ([^ with the understanding 
that results may be extended to the more general □ 

B. Equilibria and Stability Analysis by Network Reduction 

We first pursue the following question: under what condi¬ 
tions on load, network topology, admittances, and controller 
gains does the differential-algebraic closed-loop system 0 
possess a locally exponentially stable equilibrium? By exploit¬ 
ing the structure of the quadratic droop controller ([^, we will 
establish a correspondence between the equilibria of 0 and the 
solutions of a power flow equation for a reduced network. In 


this subsection we focus on generic static load models Qi{Ei) 
before addressing specific load models in Section III-C| 

For notational convenience, we first define a few useful 
quantities. We block partition the susceptance matrix and nodal 
voltage variables according to loads and inverters as 


\Eil Eli) 


E = {El,Ei), 


and define the reduced susceptance matrix 5red ^ 


( 12 ) 


by 


Bred — Ell — Bli {Bu E Kj) ^ Bjl • 


(13) 


Moreover, we define the averaging matrices Wi ^ and 

W 2 G M^x(n+m) 

^ (Bn + Ki)-^ Kj , ( 14 ) 

W2 = (Bn + Ej)-^ (-Bn Kj) . ( 15 ) 

It can be shown (Proposition |A.2| ) that Bred is invertible, and 
that Wi and W 2 are both row-stochastic matrices. Finally we 
define the open-circuit load voltages B£ G M>o ^i^ 

El ^ WiE} . ( 16 ) 


Since Wi is row-stochastic, each component of El is a 
weighted average of inverter set points B*. 

Theorem 3.1: (Reduced Power Flow Equation for 
Quadratic Droop Network) Consider the closed-loop system 
^ resulting from the quadratic droop controller ([^, along 
with the definitions The following two statements 

are equivalent: 

(i) Original Network: The voltage vector E = {El, Ei) e 

is an equilibrium point of ([8]); 

(ii) Reduced Network: The load voltage vector El 
is a solution of the reduced power flow equation 

^n = QL(EL) + [EL]Bred(EL-El), ( 17 ) 
and the inverter voltage vector is recovered via 

El = W 2 {^ij e • (18) 
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Theorem |3.1| states that, to study the existence and unique¬ 
ness of equilibria for the closed-loop system we need 
only study the reduced power flow equation ( p^ . Voltages 
at inverters may then be recovered uniquely from In 
fact, since W 2 is row-stochastic (Proposition |A.2| ), the inverter 
voltages are simply weighted averages of load voltages 
El and inverter set point voltages Ej. 

The reduced power flow equation is straightforward to 
interpret in terms of the circuit-reduction of Figure First, as 
in Remark we augment the original circuit with fictitious 
controller buses. Eliminating the inverter voltages Ej from 
the augmented network current balance ( p^ through Kron 
reduction (3T|, one obtains the input/output equivalent circuit 


II 


= J 


-^red 

/T j 


— 5red^l 

-IVi'Bred 


El 

Er 


(19) 

This reduction process is shown in Figure Left-multiplying 
the first block in ( p^ by [El] immediately yields the reduced 
power flow Hence ([HI* is exactly the reactive power 
balance equation Ql = [El\Il in the reduced network of 
Figure [Rafter reduction (cf. Equation 2.10b]). The bottom 

block of equations in ( p^ determines the fictitious controller 
current injections once the top block is solved for El. 


Proof of Theorem 3.1' (i)^(ii): Setting the left-hand side of 


■^n+m 


the closed-loop system ([8]( to zero, equilibrium points satisfy 

Ql{El) \ ^ 20 ) 

SEi]Ki{Ei - el ) ^ 

'[El] \ fB li\ (El^ 

D [Ei])\Bil Bn)\Ei^ 

Since E G solves (|2g(, we can left-multiply the lower 

block of equations in ( |20| ) by [Ej]~^ to obtain 


= Ki{Ei - E}) + BjjEj + BilEl 


( 21 ) 


and solve for Ei to obtain the inverter voltages Substitut¬ 
ing ( p^ into the first block of equations in ( [2Q| ), we calculate 


-Ql{El) = [El]BllEl 

- [EL]BLi{Bn + Ki)-^ {BilEl - KjE^ 

= [El] {Bll - BLiiBn + KiL^Bil) El 

"-V-" 

-^red 

-[El] BLi{Bn + Ki)-^KiE*i 

'■-V-" 

BredE^ 

= [El]B,,,,{El - El) , 


which is the reduced power flow equation (rf\ . 

(ii)^(i): Due to ( p^ and Proposition |A.2| (iii), we have that 
El G M>o implies that Ej G M>o hence E G M>o- 
An easy computation shows that ( p^ and ( p^ together imply 
that (El^Ei) satisfy the fixed-point equations ( [20| ), and thus 
E G is an equilibrium point. □ 

In a similar spirit, the following result states that local 
exponential stability of an equilibrium point of ^ may be 
checked by studying a reduced Jacobian matrix. 

Theorem 3.2: (Stability from Reduced Jacobian) Con¬ 
sider the closed-loop system ^ resulting from the quadratic 


droop controller •Hi. If the Jacobian of the reduced power flow 
equation given in vector notation by 

JrediEL) = ^{EL) + [EL]B,,i + [B,,LEL-El)], (22) 

is a Hurwitz matrix when evaluated at a solution El ^ M>o 
of then the equilibrium point {El, Ej) of the differential- 
algebraic system ^ is locally exponentially stable. Moreover, 
assuming that < q for each load bus i G £, a sufficient 

condition for ( |2^ to be Hurwitz is that 

Sred ^ [^l]“"[Ql(^l)]. (23) 


Proof of Theorem \3.2\ We appeal to ED Theorem 1], which 
states that local stability of the differential-algebraic system 
at the equilibrium E G may be studied by linearizing 

the differential algebraic system, eliminating the algebraic equa¬ 
tions from the system matrix, and checking that the resulting 
reduced matrix is Hurwitz. Consider the generalized eigenvalue 
problem (GEP) Jv = Xrv where v G r = blkdiag(©,r/) 
and Jacobian matrix J of 10 is given by 

J = [E]B^[BE]^D. (24) 

The diagonal matrix D G has elements Du = 

dQi{Ei)/dEi for i G £ and Du = Ki{Ei -£;*)+ KiEi for 
i ^ I. Since E is strictly positive, we left-multiply through by 
[E]~^ and formulate the previous GEP as the symmetric GEP 

Mv = X[E]-^rv, (25) 

where 


M = M^ = B + [E]-\[BE] -f D). (26) 


Partitioning the eigenvector v as v = {vl,vi), block parti¬ 
tioning M, and eliminating the top set of algebraic equations, 
we arrive at a reduced GEP Mred^/ = where 

Mred = ^rld ~ ~ ^ilM^^Mli. The matrices on both 

sides are symmetric, and in particular the matrix [Ei]~^ti 
on the right is diagonal and positive definite. The eigenvalues 
{Aijiex of this reduced GEP are therefore real, and it holds that 
< 0 for each i G X if and only if M^-ed is negative definite. 
We now show indirectly that — M^ed is positive definite, 
by combining two standard results on Schur complements. 
Through some straightforward computations, one may use (dUf, 
( p^ , ( p^ and ( [21] ) to simplify M to 


M = 


[Mil 

\Bjl 


Bli \ 

BiiEKj) ’ 


where Mil = + [^red(^L — 

El)]. Since the bottom-right block —{Bn + Ki) of —M is 
positive definite, — M will be positive definite if and only if 
the Schur complement with respect to this bottom-right block 
is also positive definite. This Schur complement is equal to 

- [El]-^ {El) - Bred - [B,ed(BL - El)], (27) 

which is exactly —[El]~^ times the Jacobian ( [22] ) of ( [TT] ). 
Since Jred is Hurwitz with nonnegative off-diagonal elements, 
— J is an M-matrix and is therefore D-stable. It follows that 
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the Schur complement ( [Z7] ) is positive definite, and hence that 
Mred is positive definite, which completes the proof of the main 
statement. For the moreover statement, proceed along similar 
arguments and consider the symmetric version 
of ( [22| ), which equals 

[El]-^^{El) + Bred + - B£)] . 

OI^L 

Solving the reduced power flow equation for 5red(^L — 
E^) and substituting into the third term above, we find that 

[EL]-^Jred = [EL]-^^iEL) + B^ed - [Bl](B^)] . 

The first term is diagonal and by assumption negative semidef- 
inite, and hence Jred(^L) will be Hurwitz if ( [2^ holds, which 
completes the proof. □ 

The sufficient condition ( [23] ) is intuitive: it states that at equi¬ 
librium, the network matrix 5red should be more susceptive 
than the equivalent load susceptances QijE];. The results of 
Theorem [3^ are implicit, in that checking stability depends on 
the undetermined solutions of the reduced power flow equation 
([Tt). The situation will become clearer in Section |III-C| when 
we apply Theorems |3.1| and |3.2| to specific load models. 

Remark 4: (Parallel Microgrids) Inwe provided an 
extensive analysis of the closed loop (jSjTfor a “parallel” or 
“star” topology, where all inverters feed a single common load. 
In particular, for a constant power load we provided a necessary 
and sufficient condition for the existence of a stable equilibrium. 
Depending on the system parameters, the microgrid can have 
zero, one, or two physically meaningful equilibria, and displays 
both saddle-node and singularity-induced bifurcations. Perhaps 
surprisingly, we show that even for a star topology, overall 
network stability is not equivalent to pair-wise stability of 
each inverter with the common load. □ 


The first technical condition in Theorem 13.31 restricts the 
impedance loads from being overly capacitive, while the second 
restricts the current loads from being overly inductive (since 
Component-wise). Both would be violated only 
under fault conditions in the microgrid, and are therefore non¬ 
conservative. As expected, in the case of open-circuit operation 
(no loading) when /shunt = ^shunt = ( [29] ) reduces to 

Ef} = E^, the open-circuit load voltage vector. 


Proof of Theorem 3.3' Substituting the 71 load model 


into the reduced power flow equation ( [Tt] ), we obtain 


— \^l\ ((-^red [^shunt])/^L (-^red-^L /shunt)) 

= [Bi](Bred + [&shunt])(BL - Ef) , (30) 


where in the second line we have factored out (5red + [^shunt]) 
and identified the second term in parentheses with (|29|). Ob¬ 


serve that El = E^^ is a solution to Since (5red + 

[^shunt]) is nonsingular and we require that El e Ef is 
the unique solution of ( [^ . It remains only to show that E'f^ G 
M>o- Since — (5red + [^shunt]) IS by assumption an M-matrix, 
the inverse (Bred + ^[shunt])”^ h^s nonpositive elements | [^ . 
By invertibility, (Bred + ^[shunt])~^ cannot have any zero rows, 
and therefore the product (Bred+ ^[shunt])“^(^redB£ -/shunt) 
with the strictly negative vector BredB£—/shunt yields a strictly 
positive vector gM”o as in ( [29] ). 

To show local stability, we apply Theorem 3.2 and check 
that the Jacobian of ( [30] ) evaluated at ( [29] ) is Hurwitz. Differ¬ 
entiating ( [^ , we calculate that 


/red — [BL](Bred + [^shunt]) + [(Bred + [^shunt])(BL “ )] , 


and therefore that 


Jred(Bf) = [Bf ](Bred + [&shunt]) , OD 

which is Hurwitz since —(Bred + [^shunt]) is by assumption 
an M-matrix and M-matrices are /J-stable. □ 


C. Equilibria and Stability Conditions for Zl, ZIP, and Dynamic 
Shunt Load Models 

In this section we leverage the general results of Section 
III-B| to study the equilibria and stability of the microgrid for 
some specific, tractable load models. To begin, the reduced 
power flow equation ( [TT] ) can be solved exactly for combina¬ 
tions of constant-impedance and constant-current loads. 

Theorem 3.3: (Stability with “ZI” Loads) Consider 
the reduced power flow equation ([Tt) for constant- 
impedance/constant-current loads 

Q^{El) = [BL][^shunt]BL + [B^J/shunt , (28) 


A more general static load model still is the ZIP model 
Qi{Ei) = 6shunt,iBf +/shunt,iBi + Qz wMch augments the ZI 
model ^2^ with an additional constant power demand (5i G M 
p5| . The model is depicted in Figure [^ In vector notation. 



^shunt,^ 
/shunt, i 


Qi 


where 6shunt ^ (resp. /shunt ^ M’^) is the vector of constant- 
impedance loads (resp. constant current loads). Assume that 

(i) -(Bred + [^shunt]) IS an M-matrix, and 

(ii) /shunt > BredB£ Component-wise. 

Then the unique solution B|^ eR% to ([Tt) is given by 

BF = (Bred + [^shunt])~^(BredB£ - /shunt) , (29) 

and the associated equilibrium point (B|^, Ef^) of the closed- 
loop system ^ is locally exponentially stable. 


Fig. 4. Pictorial representation of ZIP load model. 

the ZIP model generalizes ( [^ as 

Ql^^(Bl) = [BL][&shunt]BL + [B^J/shunt + Ql • (32) 

Unlike the ZI model ( [^ , the reduced power flow equation 
© cannot be solved analytically for ZIP loads, except for 
the special case of a parallel microgrid. Even in this special 
case, the network can have multiple equilibrium points with 

















non-trivial stability properties p9| . The following result builds 
on analysis techniques developed in pQ| and provides an 
approximate characterization of the high-voltage solution to 
(rf\ when the constant power term Ql is “small”. 

Theorem 3.4: (Stability with “ZIP” Loads) Consider the 
reduced power flow equation ( TtI with the ZIP load model ( 


let the conditions of Theorem 3.3 


hold, and let be the high- 


voltage solution of © for ZI loads as given in Theorem |3.3 
Furthermore, define the short-circuit capacity matrix 


Qsc = + [6,hunt])[^f ] 


(33) 


If IIQlII is sufficiently small, then there exists a unique high- 
voltage solution e of the reduced power flow 

equation ( p^ with ZIP loads given by 

= [Ef] (1„ - + e) , (34) 

where ||e|| = 0{\\Q~^Ql\\^), and the corresponding equilib¬ 
rium point of (U is locally exponentially stable. 

Proof: See Appendix A. ■ 

Similar to the technical conditions (i) and (ii) in Theorem [33] 
regarding constant-impedance and constant-current load com¬ 
ponents, Theorem |3.4| requires the constant-power component 
Ql of the load model to be sufficiently small. The high-voltage 
solution ( [34] ) can be thought of as a regular perturbation of the 
solution Ef^ for ZI loads, where e is a small error term. As 
IIQ^^QlII -> 0, ||e|| ^0 quadratically, and Ef}^ Ef^. 

While the ZIP load model can accurately capture the steady- 
state behavior of most aggregated loads, when considering 
dynamic stability of a power system it is often important to 
check results obtained for static load models against those 
obtained using basic dynamic load models A common dy¬ 
namic load model is the dynamic shunt model 6dyn - shunt, i = 
Qi - Efbciyn-shunt,i, or in vector notation 

^^dyn—shunt — Ql ^dyn—shunt • (3^) 


where T is a diagonal matrix of time constants. The 
model specifies a constant-impedance load model Qi{Ei) = 
bdyn-shnnt,iEf, with the shuiit susceptance 6dyn-shunt,i dynam- 
ically adjusted to achieve a constant power consumption Qi in 
steady-state. This is a common low-fidelity dynamic model for 
thermostatically controlled loads, induction motors, and loads 
behind tap-changing transformers | |35| . We restrict our attention 
to inductive loads Qi < 0, as these are the most common in 
practice, and without loss of generality we assume that < 0 
for all i G £, since if = 0 the unique steady state of ( [^ 
is 6dyn-shunt,i = 0 and the equation can be removed. 

Theorem 3.5: (Stability with Dynamic Shunt Loads) Con¬ 
sider the reduced power flow equation ( [Tt] ) with the dynamic 
shunt load model ( [^ and Qi < 0. If ||Ql|| is sufficiently 
small, then there exists a unique solution ^dyn-shunt) ^ 

^<0 of ([T7]),(|35]), given by 

i;DS = lEl] (1„ - Q-^Ql + e) , (36a) 

('dyn-shunt = [^L®]“'QL, (36b) 

where Q,, = [El]B,,d[El] and ||e|| = 0{\\Q-^QL\n and 

the corresponding equilibrium point ^dyn-shunt) 

of the system ([8]),([35]) is locally exponentially stable. 


Proof: We sketch the proof in Appendix A. ■ 

Remark 5: (From Quadratic to Conventional Droop): 

While the stability results derived above hold for the quadratic 
droop-controlled microgrid, under certain selections of con¬ 
troller gains a direct implication can be drawn from the above 
stability results to stability of the microgrid under conventional, 
linear voltage droop control. Additional information can be 
found in Section V]. □ 

We conclude this section by commenting on the utility of the 
preceding analysis and stability results. If hard design limits are 
imposed on the load voltages the expressions ( [^ or ( [34] ) 
for the unique high-voltage solutions can be used for design 
purposes to (non-uniquely) back-calculate controller gains, or 
to determine bounds on tolerable loading profiles. 

IV. Controller Performance 

Having thoroughly investigated the stability properties of 
the closed-loop system in Section [I^ we now turn to 
questions of optimality and controller performance. Here we 
are interested in inverse-optimality of the resulting equilibrium 
point of the system in an optimization sense, along with 
how well the resulting equilibrium achieves the power sharing 
objective. We do not addresses system-theoretic optimality or 
performance of the controller. 


A. Optimality of Quadratic Droop Control 


While in Section III we emphasized the circuit-theoretic 
interpretation of the quadratic droop controller ([^, it is also 
relevant to ask whether the resulting equilibrium point of 
the system is inverse-optimal with respect to any criteria. 
For simplicity of exposition, we restrict our attention in this 
subsection to constant-impedance load models of the form 


Qi{^i) — ^shunt,i-^i 


(37) 


Analogous extensions to ZIP load models are possible using 
energy function theory | [42| , at the cost notational complex¬ 
ity and less insight. To begin, note that for nodal voltages 
{El, Ei) e the total reactive power Qioss absorbed by 

the inductive transmission lines of the microgrid is given by 

Q\oss{El^ Ej) = Bij{Ei — Ej) . 

Similarly, the total reactive power consumed by the loads is 

Qload(-^L) — ^ ^ . ^shunt,i-^2 • 

Efficient operation of the microgrid stipulates that we minimize 
reactive power losses and consumption, since transmission and 
consumption of reactive power contributes to thermal losses 
Simultaneously however, we have the practical requirement that 
inverter voltages Ei must remain close to their rated values 
E’* >0. We encode this requirement in the quadratic cost 

Cvoit {Ei) = -y^ Ki{Ei-E*)‘^, (38) 


^For an inverse-optimality analysis that relates real power generation to 
frequency droop control, see 0. 
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where K,i < 0 are cost coefficients. Consider now the combined 
optimization problem 


minimize C = Qioss{E) + Qioad(-E’L) + C'voit(-E’7), (39) 


where we attempt to enforce a trade-off between minimizing 
reactive power dissipation and mininfizing voltage deviations. 
The next result shows that an appropriately designed quadratic 
droop controller •Hi is a decentralized algorithm for solving the 
optimization problem ( [39] ). Conversely, any quadratic droop 
controller solves an optimization problem of the form ( [^ for 
appropriate coefficients Ki. 

Proposition 4.1: (Optimality and Quadratic Droop) Con¬ 
sider the closed-loop microgrid system with constant- 
impedance loads and controller gains {Ki}i^x^ and the 
optimization problem ([3^ with cost coefficients 
Assume as in Theorem |3.3| that —{Bred + [^shunt]) is an M- 
matrix. If the parameters are selected such that Ki = Ki 
for each i e I, then the unique locally exponentially stable 
equilibrium {Ef^.Ef) G of the closed-loop system ^ 

is equal to the unique minimizer of the optimization problem 
( [39| ), and both are given by 

El = (Bred + [&shunt])“'Bred^I € , (40) 


with Ef e 


Proof of Proposition 4.1 


>0 recovered by inserting ( [40] ) into 

That {Ef, Ef) is the unique locally 


exponentially stable equilibrium of ^ follows immediately 
from Theorem 33 by setting /shunt = We now relate the 
critical points of the optimization problem ( [^ to the equilib¬ 
rium equations for the closed-loop system. In vector notation, 
we have that Qioss — E BE, Qioad ~ [^shuntj-E'L? 
and (3voit = —{Ej — Ej)^Kj{Ei — Ej), where E = 
{El, ..., En-\-m)’ The total cost C may therefore be written as 
the quadratic form C = x~^Bx, where x = {El, Ej, Ej) and 


/ Bll + [^shunt] Eli 0 \ 

S = - Bil Bn + Ki -Kj 

\ 0 -Ki Ki ) 


The first order optimality conditions VC{E) = ©n+m yield 

fB ll E[hshnnt] BlI A (El\ ©n A /..X 

Bil Bii^Ki)\Ei)-\KiE^j)' 

Solving the second block of equations in ( |4T] ) for Ei yields 
the quadratic droop inverter voltages Substituting this 
into the first block of equations in ( [4T] ), simplifying, and left- 
multiplying by [El] yields 


©n = [EL][bshunt]EL + [/^L]^red {El “ Ef) , (42) 

where Bred and Ef are as in ([T^-([T^. The result is 
exactly the reduced power flow equation ( [30| ) with /shunt = 
©n, which shows the desired correspondance. It remains only 
to show that this unique critical point is a minimizer. The 
(negative of the) Hessian matrix of C is given by the matrix of 
coefficients in Since the bottom-right block of this matrix 
is negative definite, it follows by Schur complements that the 
Hessian is positive definite if and only if —(Bred + [^shunt]) 
is positive definite, which by assumption holds. □ 


Note that depending on the heterogeneity of the voltage set 
points B*, minimizing Qioss + Qioad may or may not conflict 
with minimizing Cvoit- The quadratic droop controller is a 
primal algorithm for the optimization problem ( [3^ , and can 
therefore be interpreted as striking an optimal balance between 
maintaining a uniform voltage profile and minimizing total 
reactive power dissipation (cf. | [T9| , | [43| , where similar trade¬ 
offs are studied). The result can also be interpreted as an 
application of Maxwell’s minimum heat theorem to controller 
design Proposition 3.6]. 

Remark 6: (Generalized Objective Functions) To obtain 
the generalized feedback controller ( [TT] ), one may replace 
C'voit(£^/) given in @ by 

C.oii{Ei) = V, . ^ "^EiEj + pijEiE*. 

For a general strictly convex and differentiable cost C^o\t{Ei), 
analogous methods may be used to arrive at feedback con¬ 
trollers of the form Ui = E^ • See |j^ Chapter 9], 

Remark 4] pT| , (2^ for more on nonlinear droop curves 
derived from optimization methods, including dead bands and 
saturation. □ 


B. Power Sharing: General Case and Asymptotic Limits 


We now study the following question: under an incremental 
change in load power demands Ql Ql + how do 
the reactive power injections at inverters change? This is the 
question of power sharing. For simplicity of exposition, we 
make two standing assumptions in this subsection. First, we 
consider only the case of constant power loads Ql{El) = Ql- 
Second, we consider the case of uniform inverter voltage set 
points, such that = Bat • for some fixed nominal voltage 
Bat > 0. Both assumptions can be easily relaxed to yield more 
general (but much more cumbersome) formulae. Under these 
assumptions, the results of Theorem |3.4| simplify to yield the 
approximate load voltages 


-El - -Biv ( In - -f^B\QL 


B^ 


(43) 


with Bred as defined in Theorem 3.1 The following definition 
is useful for interpreting our result to follow. 


Definition 1: (Differential Effective Reactance) Let E^^ 

be the open-circuit voltage at load i ^ C. The differential 
effective reactance between loads G £ is the voltage 
difference Ei — Ef^ measured at load i when a unit current is 
extracted at load j with all other current injections zero. 

In other words, is the proportionality coefficient from 
current injections at load j to voltage deviations at load i. 
For branches {i,j} G £, we also let Xij denote the direct 
reactance between buses i and j, and set Xij = +00 if there 
is no branch. The main result of this section characterizes how 
inverter reactive power injections behave under the quadratic 
droop controller ([^ as a function of the susceptance matrix B, 
load demands Ql, and controller gains Ki. 


Theorem 4.2: (Power Sharing) Consider the closed-loop 
system © with the approximate load voltages ( [4^ . Then the 
inverter reactive power injections Q/ G are given by 


Ql = KiiBii + Ki)-^BilB;2Ql ■ (44) 
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Moreover, the following special cases hold: 

1) Distance-Based Power Sharing: In the high-gain limit 
where Ki ^ —oo for each inverter i G X, it holds that 


Qi 



i el. 


(45) 


That is, loads are preferentially supplied with power by 
sources which are electrically nearby; 

2) Proportional Power Sharing: In the low-gain limit 
where Ki ^ 0 ~ for each i el, it holds that 

~ -^Qtotai: i e I, (46) 

where Qtotai = Qi losid. That is, sources 

supply power based only on their relative controller gains. 


In the high gain limit ( [45] ), each inverter behaves as a 
stiff (constant) voltage source, and power is routed based on 
electrical distances. In particular, ( [45] ) adds up the parallel 
transfer paths from all loads j through all intermediate loads 
k adjacent to inverter i. In the low-gain limit of ( |46| ), inverters 
react to changes in total load by supplying reactive power in 
proportion to their controller gain Ki, which can be selected 
as being proportional to the generation capacity of the unit. 
Moreover, observe that ( |4^ does not depend on the network 
topology. The general formula ([g can be interpreted as 
interpolating between these two regimes. 

Proof of Theorem \4.2\ From the closed-loop system 
we observe that in steady-state, the inverter reactive power 
injections are given by 


Qi = Ki[Ei]{Ei — , (47) 


in which the only variables are the inverter voltages Ej. 
However, Theorem 3T] yielded an expression for Ej in terms 
of El, namely ([18). Substituting the load voltages into 
and working through some algebra, one finds that Ej = 
En • W^^tm AiQl, where W 3 is row-stochastic, and hence 


El = Eulm + AiQl , (48) 


where 

Ai ^ -f(Bn + . (49) 

rjN 

Much like in the proof of Theorem |3.4| we linearize ( [47] ) around 
the open-circuit operating condition {Ej = Entm^Qi = ©m)- 
Performing this linearization, we obtain 

Qi = EnKi{Ei — E]s['^m ), 

which after substituting ( |48] ) and ( [4^ yields the desired result 
(|44|). To show the first statement, for 5 > 0 define 

Bved,E = Bll — Bli{Bii -h eKi)~^BiL (50a) 
S, = sKiiBii + sKi)-^BilB-1^ . (50b) 

Note that for 5 = 1, we recover the previously derived formula 
( |44| ) as Q/ = S'iQl- We calculate that lim^^oo ^red,£ = Bll, 
and that lim^^oo ^Ki{Bii + eKi)~^ = /^. It follows that 

lim = BjlBj^I . 


Applying Lemma [A3| then leads to the first statement. The 
proof of the second statement is more delicate, since for 5 = 0, 
5red,o = Bll—BliBj}Bli is singular with one-dimensional 
kernel spanned by tm Lemma II. 1]. Nonetheless, 5“^ ^ 
may be expanded in a L^rent series around the simple pole 
at 5 = 0 |[^. Moreover, since lim^^o = 

eKjBj^, it is clear that the only term in the Laurent series of 
^red e ^Eing rise to a non-zero contribution to lim^^o is the 
term proportional to e~^. Some straightforward computations 
using | |44| Theorem 1.1] give that the required term in the 
Laurent series is and hence that 


limS'e = lim sKjBjIBil ■ - J , Inxn , 


1 




KiBj^BiLt 


mxn • 


Since —BjjBjl is 
Bii BiL^nxn — 
things out in components using Lemma [A3 


row-stochastic |[31j Lemma 11.1], 
mxn, and the result follows by writing 

□ 


As a final technical remark regarding Theorem 4.2 we 
emphasize that the proportional power sharing result 
represents a limit under the linear approximation (|43. In 
this low-gain limit, 5“^ in becomes very large and 
for any fixed Ql the second term in the linearization ( [^ 
dominates, leading to large voltage deviations. This violates 
the premise under which the linearization was derived, and 
may lead to the loss of the closed-loop system’s equilibrium 
point or to violation of the decoupling assumption. Thus, the 
extent to which the proportional power sharing result may be 
implemented in practice is limited by the size of the loads to be 
serviced, the stiffness of the grid, and the stability bottleneck 
associated with low load voltages. We explore this issue further 
via simulations in Section m 

Nonetheless, an interesting observation can be made by 


comparing Theorem 4.2 with the optimization problem 
Namely, large controller gains lead to accurate voltage regu¬ 
lation and distance-based reactive power sharing, while small 
controller gains minimize reactive power dissipation and yield 
proportional power sharing (but with risk of instability). 


V. Simulations 

Through simulations we now demonstrate the robustness 
of our results to unmodeled resistive losses, coupled ac¬ 
tive/reactive power fiows, and frequency dynamics. We consider 
the islanded microgrid of Figure [^b) consisting of five loads 
and three inverters. Loads are modeled using the dynamic 
shunt model for reactive power demands, along with 
the analogous dynamic conductance model for active power 
demands. The consumption set points in these load models 
are taken as disturbance inputs in what follows, in order to 
subject the system to time-varying loads. We model additional 
high-frequency variation on the load consumption set points 
as uncorrelated white noise, with a standard deviation of 
15% of the nominal set points for both active and reactive 
power. Branches have a mixed resistive/inductive character 
with non-uniform R/X ratios of between 0.3 and 1.1. Inverter 
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frequencies are controlled by the frequency droop controller 

cj- = cj* - rriiPe^i , i G X, 

where uji is the inverter output frequency, cj* = 27r • 60Hz, and 
mi > 0 is the frequency droop coefficient. System parameters 
are based on those from the experiments in Q, and are omitted 
for brevity. The quadratic droop gains are set based on the 
conventional droop gains from |7| as which is 

a relatively stiff tuning leading to strong voltage regulation. In 
particular, Ki = K 2 = meaning that inverters 1 and 2 
are rated equally, with inverter 3 being rated for half as much 
power. 

The system is under nominal loading up to t = 2s, at which 
point loads {1, 3, 5} are subject to sinusoidal variations with a 
50% magnitude and a 0.5s period for both active and reactive 
power (at constant power factor). At t = 4s sinusoidal variation 
stops and active and reactive consumption set points at load 
buses {2,4} are doubled, before returning to nominal loading 
at t = 6s. Figure shows the resulting inverter reactive power 
injections!^ With this tuning, the controller is able to maintain 
stability and track the load demands, but proportional power 
sharing is poor in accordance with Theorem |4.2| (the curves 
for inverters 1 and 2 should overlap, but do not). Figure 
shows the resulting voltage magnitudes of both inverter and 
load buses. The inverter voltages stay well regulated, which 
prevents load buses voltages from falling even lower. 


Inverter Reactive Power Injections 



Fig. 5. Evolution of inverter reactive power injections under quadratic droop 
control for large feedback gains K^. Vertical axis normalized by base power 
of 1400 VA. 

We now test the prediction from Theorem |4.2| that making 
Ki small will enforce proportional power sharing, but also make 
the network more prone to instability. We reduce the feedback 
gains Ki substantially, to 5% of their previous values, and 
repeat the above experiment. The resulting traces are shown in 
Figures and As expected, early in the experiment reactive 
power sharing is enforced, with inverters 1 and 2 sharing power 
equally and providing twice the power of inverter 3. As seen in 
Figure [Sj voltage magnitudes are significantly lower than they 

§The corresponding traces under the conventional voltage droop controller 
are extremely similar, and are therefore omitted for clarity. 


Nodal Voltage Magnitudes 



Fig. 6. Evolution of voltage magnitudes under quadratic droop control for 
large feedback gains K^. Vertical axis normalized by base voltage of 230V. 


were under higher controller gains. When the load is doubled 
at t = 4s a voltage collapse occurs and the system is unable 
to recover. This confirms that our theoretical results remain 
useful in a more general model setting. 


Inverter Reactive Power Injections 



Eig. 7. Evolution of inverter reactive power injections under quadratic droop 
control for small feedback gains Ki . Vertical axis normalized by base power 
of 1400 VA. 


VI. Conclusions 

In this work we have provided a detailed analysis of the 
quadratic droop controller, a tractable modification of conven¬ 
tional droop control for application in islanded microgrids. 
The special structure of quadratic droop control allows for 
rigorous circuit-theoretic analysis methods to be applied to the 
closed-loop system. A detailed analysis of system equilibria 
and stability was provided for several load models, along with 
an optimization interpretation of the resulting decentralized 
controller. The power sharing characteristics of the controller 
were analyzed, showing that the controller interpolates between 
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Nodal Voltage Magnitudes 



Fig. 8. Evolution of voltage magnitudes under quadratic droop control for 
small feedback gains Ki. Vertical axis normalized by base voltage of 230V. 


proportional power sharing (in the low-gain regime) and power 
sharing based on electrical distance (in the high-gain regime). 
Our analysis results provide easily verifiable certificates for 
system stability, and the proof techniques developed should 
prove quite useful for related microgrid analysis and control 
problems. 

The approach and techniques reported in this paper should 
serve as a guidepost for further work on the following key 
technical problems. While our results apply to microgrids 
with uniform R/X ratios under a decoupling assumption, the 
coordinate transformation used herein to treat uniform R/X 
ratio microgrids (see Remark has the side-effect of mixing 
active and reactive power, making true active/reactive power 
sharing more subtle to implement. Relaxing these assumptions 
may prove possible using recently developed analysis tools 
such as coupled power flow linearizations. The design of a 
simple, decentralized, and provably stable droop-like controller 
for voltage stabilization and power sharing in microgrids 
with non-uniform R/X ratios remains an open problem. 
Additional important practical directions are an analysis of 
interactions between inner voltage/current control loops and 
droop control, along with a voltage stability analysis for mixed 
inverter/generator networks. 

While here we have addressed the problem of primary 
control, the authors view secondary control for islanded mi¬ 
crogrids 0 as an area requiring new problem formulations 
and increased theoretical attention. In particular, it is not clear 
to what extent standard secondary control formulations such 
as voltage regulation and reactive power sharing (adopted 
from high-voltage systems and power supply backup systems, 
respectively) are important in sub-distribution sized microgrids. 
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(i) Bij = Bji > 0 for i 7 ^ j, with strict inequality if and 
only if {i,j} e £; 

(ii) B is negative semidefinite with a simple eigenvalue at 
zero corresponding to the eigenvector ln-\-m ; 

(iii) all principal submatrices of B are negative definite. 


Proposition A. 2: (Properties of Reduced Quantities) The 
following statements hold: 

(i) —^red is an M-matrix; 

(ii) Wi and W 2 are row-stochastic; 

(iii) E^l > component-wise. 


Proof of Proposition A. 2“ (i): The first fact follows from the 


closure of the set of symmetric positive definite M-matrices 
under the Schur complement, as 5red is the Schur complement 
of 5 + blkdiag(©, Kj) with respect to the Bn block fsT] 
Lemma II. 1]. Similar to (U Lemma II.l], it can be verified 
that Wilm = In and IL 2 In+m = Im- The proof of (iii) then 
follows from (ii). □ 


Proof of Theorem 3.4' With the ZIP load model the 


reduced power fiow equation ([30|) is modified to 


Ql = - [El\{Byq^ + [^shunt])(^L - • (51) 

Consider now the invertible change of coordinates defined by 

EL = [Ef]{t,,-Q-,^QL + e), (52) 

where e G is our new variable. Substituting this into the 
power fiow (Hg, some basic algebra leads to 


^(e, Ql) = QscE + [qL - e]Qsc(9L - e) = ©n , (53) 

where we have abbreviated ql = Q^IQl- First, note 
that /i(©n,©n) = ©n. Next, an easy calculation gives that 
H(©^, ©^) = Qsc, which is full rank. It follows from the Im¬ 
plicit Function Theorem that there exist open sets f/q, kq C 
(each containing the origin) and a continuously differentiable 
function H : Uq ^ Vq such that h{H{qL)^qL) = fof 
all qp G Uq. This in turn shows that for 119^1,11 = ||Q“^Ql|| 
sufficiently small, the expression ( [52| ) solves the reduced power 
flow equation ( [ST] ). To obtain the quadratic bound on ||e||, 
rearrange ( [5^ and take the norm of both sides to obtain 

IIQscell = \\[qL - e]Qsc{qL - e)ll • (54) 

Since Qsc is invertible, the left-hand side of •HI is lower- 
bounded by <T||e||, where a > 0. Moreover, a series of 
elementary inequalities shows that the right-hand side of ( [54| ) 
is upper bounded by (3\\qL — e|p < /5(||gL|| + IklD^^ where 
P > a depends only on Qsc- Combining these bounds, we find 
that e satisfies 


Appendix A 

Supporting Lemmas and Proofs 


||e||<^(||9L|| + ||e||)L 
a 


(55) 


Lemma A.l: (Properties of Susceptance Matrix) The sus- 
ceptance matrix B G satisfies 


Since WqpW is sufficiently small, in particular it holds that 
dallq'LlI//^ < 1. In this case, some simple algebra shows that 
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inequality ( [55] ) further implies that 

IWI < ^ - Ikill - 

= d\\qa^ = 0{\\Q-^Qd^), 

where we have used the fact that for a G [0,1], Vl — a > 
1 — a/2 — a^/2. Exponential stability follows immediately by 
continuity from the case of ZI loads. □ 


Proof of Theorem |5.5[ We provide only a sketch of the 

is equivalent to 


proof. In steady state, the DS model 
a constant-power model, and hence the approxi mate steady- 
state load voltages are given as in Theorem 3.4 with ZI 
components 6shunt = and /shunt = yielding the stated 
equilibrium. Stability can no longer be shown by applying 
Theorem |3.2[ since Theorem |3.2| was proven in the absence of 
load dynamics. A tedious but straightforward calculation shows 
that applying the proof methodologies of Theorems |3.2| and 
|3.4| to the extended dynamics ([8]),([^ yields local exponential 
stability. □ 
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Lemma A. 3: (Electrical Distances) Let be the differ¬ 
ential effective reactance between loads k^j G C, and for 
{i^k} G E let Xik = l/Bik be the direct reactance, with 
Xik = +OC if {i, k} ^ E. Then 


-{BjlB 


Tl) ij 




Proof of Lemma |A./[ From the current balance equations 
—]Il = BllElA-BliEj, observe by setting II = ©n that the 
open-circuit load voltages are given by = —BflBLiEj. 
Then the current balance can be rewritten sls El — EJf = 
—jBfllL, and it follows that X^^ = —Bf\. The desired 
formula is then immediate by noting that for i G £, j G X, 
and {i, j} G E, ll{BLi)ij = Xij by definition. □ 
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